Spin- Wave Modes and Their Intense Excitation Effects in Skyrmion Crystals 
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We theoretically study spin-wave modes and their intense excitations activated by microwave 
magnetic fields in the skyrmion-crystal phase of insulating magnets by numerically analyzing a 
two-dimensional spin model using the Landau-Lifshitz-Gilbert equation. Two peaks of spin-wave 
resonances with frequencies of ~1 GHz are found for in-plane a.c. magnetic field where distribution 
of the out-of-plane spin components circulates around each skyrmion core. Directions of the circula- 
tions are opposite between these two modes, and hence the spectra exhibit salient dependence on the 
circular polarization of irradiating microwave. A breathing-type mode is also found for out-of-plane 
a.c. magnetic field. By intensively exciting these collective modes, melting of the skyrmion crystal 
accompanied by a red shift of the resonant frequency is achieved within nano seconds. 

PACS numbers: 76.50.+g, 75.70.Ak, 75.10.Hk, 75.78.-n 



Competing interactions in magnets often cause non- 
trivial spin textures such as ferromagnetic domains and 
magnetic bubbles, which have attracted a great deal of 
interest from the viewpoints of both fundamental science 
and technical applications in the field of spintronics [HQ- 
In particular, response dynamics of such magnetic struc- 
tures under external fields is an issue of vital importance 
because its understanding is crucial for their manipula- 
tions. 

Skyrmion, a nontrivial swirling spin structure carry- 
ing a topological quantum number, is one of the inter- 
esting examples of such spin textures. It was originally 
proposed by Skyrme to account for baryons in nuclear 
physics in 1960's as a quasiparticlc excitation with spins 
pointing in all directions to wrap a sphere 0, and 
was recently realized experimentally in two-dimensional 
condensed matter systems, e.g., quantum Hall ferromag- 
nets 0, [(§ , ferromagnetic monolayers , and doped lay- 
ered antiferromagnets Q- 

Formation of skyrmion crystal (SkX) was theoreti- 
cally predicted in Dzyaloshinskii-Moriya (DM) ferromag- 
nets without inversion symmetry 0, [l(|, and was in- 
deed observed in the A-phase of metallic chiral magnets 
MnSi [ll, 13 [ and Fei-^Co^Si [l2j by neutron-scattering 



experiments as a triangular lattice of skyrmions with 
spins antiparallel to the applied magnetic field at the 
skyrmion centers and parallel at their peripheries. A re- 
cent Monte-Carlo study found a greater stability of the 
SkX phase in thin films [l4[. This prediction was con- 
firmed by the real-space observation of skyrmion trian- 
gular lattice in Feo.5Coo.5Si thin films using the Lorentz 
force microscopy in a wide temperature and magnetic- 
field range 15 1. 

Typically the skyrmion is 10-100 nm in size, which is 
determined by the ratio of DM interaction and exchange 
coupling and is much smaller than magnetic bubbles. 
Moreover, recent experiments found that the skyrmion is 
stable even near/above the room temperature fltj ]. and 
can be manipulated by much lower electric currents than 



ferromagnetic domain walls [171 lla|- These properties, 
i.e., small size, high operational temperature, and low 
threshold field, are advantageous for technical applica- 
tion to high-density data storage devices. Therefore, un- 
derstanding of the dynamics of skyrmions and SkX under 
external fields is an important issue fl9j |. 

In this Letter, we theoretically study collective spin dy- 
namics in the SkX phase of insulating fcrromagncts with 
DM interaction by numerical simulations of the Landau- 
Lifshitz-Gilbert (LLG) equation under time-dependent 
a.c. magnetic fields. We find a couple of spin-wave reso- 
nances with frequencies ~ 1 GHz for in-plane a.c. mag- 
netic field where the out-of-plane spin components rotate 
around each skyrmion core. The directions of these rota- 
tions are opposite between the higher-lying and lower- 
lying modes, and their spectra show strong circular- 
polarization dependence. A breathing-type mode is also 
found for out-of-plane a.c. magnetic field. Furthermore, 
we study intense excitation effects of these collective 
modes, and find a red shift of the resonant frequency 
and melting of the SkX within nano seconds. These find- 
ings will lead to a fast manipulation of skyrmions in nano 
scale using spin- wave resonances. 

We start with a classical Heisenberg model on a two- 
dimensional square lattice (l4j |. which contains nearest- 
neighbor ferromagnetic exchange, Zeeman coupling, and 
DM interaction as 25|, 



H 



-J Y, Si-Sj-iH + H'W-YSi 



<i,]> 



Dj2(Si x s i+x -x + SiX S i+ y ■ y), (1) 



where H=(0, 0, H z ) is a constant external magnetic 
field normal to the plane, and H'(t) is an applied time- 
dependent magnetic field. The norm of the spin vector 
is set to be unity. We adopt J=l as the energy unit and 
take £)=0.09. The spin turn angle 9 in the helical struc- 
ture is determined by the ratio D/J&s tan#=£)/(\/2J), 
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FIG. 1: (color), (a) Phase diagram of the Hamiltonian |T|) 
at T=0 where HL, SkX, and FM denote helical, skyrmion- 
crystal, and ferromagnetic phases, respectively, (b) Spin con- 
figuration of the SkX phase with a color map of the spin z- 
axis components S z i at H z =3.75xl0 _3 . Spin vectors at sites 
(i x , i y ) projected onto the xy plane are shown by arrows for 
mod(ia;, 6)=mod(i !/ , 6)=0. (c) One skyrmion is magnified 
with a color map of the scalar spin chiralities d. 



TABLE I: Unit conversion table when J— OA meV. 



magnetic field H 


lxl(T 3 J 


-3.4 mT 


frequency u 


0.01J 


~1 GHz 


time t 


toooj- 1 


-10 nsec 



which is derived from a saddle point equation of the en- 
ergy as a function of 9. Our parameter set gives 0=3.64° 
or the periodicity of ~99 sites, which corresponds to the 
skyrmion diameter of ~50 nm if we consider a typical 
lattice parameter of 5 A. 

We study collective spin excitations of this model by 
numerically solving the LLG equation using the fourth- 
order Rungc-Kutta method. The equation is given by 



dSj, 
dt 



1 



Si x H 



off 



s 



Si x {Si x Hf) 



(2) 

where ac is the dimcnsionlcss Gilbert-damping coeffi- 
cient. We derive a local effective field Hf acting on 
the ith spin Si from the Hamiltonian H as Hf = 
—dH/dSi. All the calculations are performed for sys- 
tems with Af=288x288 sites under the periodic boundary 
condition. We fix aQ=0.04 for simulations of the spectra 
shown in Fig. [2j while o:g=0.004 for others. 

We first study phase diagram of the model (JT|) at T=0 
as a function of H z . Starting with spin configurations 
obtained in the Monte-Carlo thermalization at low T, 
we further relax them by sufficient time evolution in the 
LLG equation, and compare their energies. As shown 
in Fig. QJa), helical (HL), SkX, and ferromagnetic (FM) 
phases appear successively as H z increases where criti- 
cal fields are 7J 2 =1.875x 10 -3 and H z =6.3x 10 -3 , respec- 
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FIG. 2: (color online). Imaginary parts of (a) in- plane and (b) 
out-of-plane dynamical susceptibilities, Imx(w), in the SkX 
phase for several values of H z . Insets show resonant frequen- 
cies wh as functions of H z . 



tively. Here H z =lxl0~ 3 corresponds to ^3.4 mT if we 
adopt a typical value of J=0.4 meV and 5=1 spins (see 
also Table |T|). In Fig.QJb), we display spin configuration 
of the SkX phase where the in-plane components of the 
spin vectors at sites (i x , i y ) are described by arrows when 
mod^x, 6)=mod(ij,, 6)=0. Here distribution of the spin 
z-axis components, S Z i, is shown by a color map. One 
skyrmion is magnified in Fig. [ljc) with a color map of 
the local scalar spin chiralities given by, 

C'i = Si ■ (Si + x X Si + y) + Si ■ (Si — X X Si — y). (3) 

The finite spin chirality is a source of the topological Hall 
effect observed in experiments 21 - 24 1 . 

We then study the microwave-absorption spectra due 
to spin-wave resonances in the SkX phase. We trace spin 
dynamics after applying a (5-function pulse of magnetic 
field at t=0, which is given by H'(t)=S(t)H u . The ab- 
sorption spectrum or the imaginary part of the dynami- 
cal susceptibility, Imx(w), is calculated from the Fourier 
transformation of magnetization m(t)=(l / TV) J2i Si(t). 

In Fig.[5Ja), we show calculated spectra for several val- 
ues of H z when H^ is parallel to the xy plane. We find 
two resonance peaks in the spectra, and both of their fre- 
quencies increase as H z increases as shown in the inset 
of Fig. [21a). Note that w=0.01 corresponds to ~1 GHz 
for ,7=0.4 meV (=96.7 GHz). Thus these spin- wave res- 
onances are located in the frequency range 500 MHz-1.2 
GHz or in the microwave regime. On the other hand, 
the calculated spectra for H^ parallel to the z axis are 
shown in Fig. EJ^b) , which have only one resonance peak. 
The resonant frequency ojr decreases as H z increases as 
shown in the inset. Again these resonances are located 
in the microwave frequency regime. 

To identify each spin- wave mode, we trace the spin dy- 
namics by applying a stationary oscillating magnetic field 
with resonant frequency cjr. We first study the modes 
activated by the in-plane a.c. magnetic field by set- 
ting H'(t)=(0, H"smijj K t, 0) with i^=0.5xl0~ 3 . The 
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FIG. 3: (color). Spin dynamics of each collective mode in the SkX phase calculated at ff 2 =3.75xl0 -3 . Spins at sites iy) 
are represented by arrows when mod(ia;, 6)=mod(i H , 6)=0 with color maps of the S Z i components in the left panels, while in the 
right panels, distributions of the local spin chiralities G% are displayed. Temporal waveforms of the applied a.c. magnetic fields, 
Hy sintJRt and H% smuiRt, are shown in the uppermost figures where inverted triangles indicate times at which we observe the 
spin configurations shown here, (a) [(b)] Lower-energy [Higher-energy] rotational mode with ojr=6.12x10 _3 (o;r=1.135x10 -2 ) 
activated by the in-plane a.c. magnetic field. Distributions of the S Z i components and the spin chiralities d circulate around the 
skyrmion core in a counterclockwise (clockwise) fashion, (c) Breathing mode with u;r=7.76x 10 -3 activated by the out-of-plane 
a.c. magnetic field. 



frequency ur is fixed at wr=6.12x 10~ 3 for the lower- 
energy mode, while at wr=1.135x 10~ 2 for the higher- 
energy mode. We find that for all of the modes, all the 
skyrmions show uniformly the same motion so that we 
focus on one skyrmion hereafter. In Figs. EJa) and (b), 
we display calculated time evolutions of the spins. The 
spins at sites (i x , i v ) are represented by arrows when 
rao&(i x , 6)=mod(ij,, 6)=0 together with distributions of 
the S Z i components in the left panels, while those of the 
spin chiralities d in the right panels. Interestingly the 
area of larger S z i or that of larger \d\ circulates around 
each skyrmion core even though the applied a.c. field 
H'(t) is linearly polarized in the y direction. We find 
that directions of their rotations are opposite, i.e., coun- 
terclockwise (CCW) with respect to the magnetic field 
H\\z for the lower- lying mode while clockwise (CW) for 



the higher-lying mode. These directions are independent 
of the sign of DM constant D or winding direction of 
the spins. Instead they are determined by a sign of the 
applied field or by the spin orientation at the skyrmion 
core. 

Because of these habits, the spin- wave excitations acti- 
vated by the in-plane a.c. magnetic field strongly depend 
on the circular polarization of the irradiating microwave. 
In Fig.lH we show calculated time evolutions of the mag- 
netization parallel to the y axis, m y (t)={\ / N) J2i S y i(t), 
when we irradiate linearly-polarized, left-handed circu- 
larly polarized (LHP), and right-handed circularly po- 
larized (RHP) in-plane microwaves with resonant fre- 
quency wr=6.12x 10~ 3 , which corresponds to the lower- 
lying mode at ff z =3.75x 10 -3 . More concretely, we apply 
a time-dependent magnetic field H'(t)=\H' x (t), H'(t), 0] 
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FIG. 4: (color online). Calculated time evolutions of mag- 
netization (||y), m y (t)=(l/N) J]V S y i(t), in the SkX phase 
at 7^2=3.75 xlO -3 under linearly-polarized, left-handed circu- 
larly polarized (LHP), and right-handed circularly polarized 
(RHP) in-pane a.c. magnetic fields with resonant frequency 
cjr=6.12x 10 -3 corresponding to the lower-lying mode. 



(b) i=4000 m "•. (c) f=7200 ■ 




FIG. 5: (color online). Melting of the SkX within nano sec- 
onds under irradiating LHP microwave, which excites the ro- 
tational spin- wave modes intensively (see text). Color maps of 
the S z i components are displayed at (a) t=0, (b) i=4000, and 
(c) t=7200. Figures magnify a partial area with 220x220 sites 
for clarity, while the calculations are done for 288x288 sites 
with the periodic boundary condition. Temporal waveform 
of the microwave is also shown where times corresponding to 
figures (a), (b), and (c) are indicated by inverted triangles. 



where H' x {t)=aH xy coscjRi and H' (t)=H£ sinwRi with 
a=0 for the linearly-polarized microwave and a=l ( — 1) 
for the LHP (RHP) microwave. In the LHP (RHP) mi- 
crowave, its magnetic-field component rotate in a CCW 
(CW) way. Here we fix H% y =0.5x 1CT 3 . We find that 
irradiation of the LHP microwave significantly enhances 
the magnetization oscillation as compared to the linearly- 
polarized microwave, whereas the RHP microwave can- 
not activate collective spin oscillations. 

Next we discuss a spin-wave mode activated by the 
out-of-plane a.c. magnetic field. We again trace 
spin dynamics by applying H'(t)=(0, 0, H" sin wr£) 
with w R =7.76xlCr 3 and H^=0.5x 10" 3 . We observe a 
breathing mode where the area of each skyrmion extends 
and shrinks dynamically as shown in Fig. |3]Jc) . 

We finally study effects of the intense spin-wave exci- 
tation. We apply in-plane LHP (a—1) and RHP (a=— 1) 



microwaves of H' x (t)~aH xy coswi and H' y (t)=H xy sinwi 
to the SkX phase at H z ~6.3x 10~ 3 . The system is lo- 
cated on the phase boundary between the SkX and FM 
phases. Here we take H xy =0.5x 10~ 3 , which corresponds 
to ~1.7 mT when J=0.4 meV and 3=1. The frequency 
u is fixed at 7.4xl0 -3 . This value is nearly equal to 
the resonant frequency wr=7.8x 10 -3 of the lower-energy 
mode, but slightly deviates from it in reality. Because 
the intense spin-wave excitations necessarily change the 
spin structure from its equilibrium configuration, and 
it results in red shifts of the resonant frequencies, we 
chose lo slightly smaller than ojr of the nearly equilib- 
rium case in advance. In fact, the red shift can be seen 
in Fig. |U The magnetization dynamics under the LHP 
microwave becomes slow as compared to that under the 
linearly-polarized microwave when the oscillation ampli- 
tude becomes larger. One can easily notice this fact 
from different maximum points between these two oscil- 
lations. Indeed the oscillation frequency in Fig. blunder 
the LHP microwave is w^6.1xl0 -3 for 0<t<2000, while 
w~5.7xl(r 3 for 3000<t<5000. 

In Figs. [5ja) and (b), we show snap-shots of the spin 
configurations at several times under the irradiating LHP 
microwave. We observe melting of the SkX due to the in- 
tensively excited rotational spin-wave modes. The melt- 
ing occurs within £^5000-6000. Here £=1000 corresponds 
to ~10 nsec when J=0.4 meV. Thus the melting occurs 
within 50-60 nsec. We also find that the SkX melting is 
difficult to achieve either by the RHP microwave or even 
by the LHP microwave if its frequency is off- resonant. 
Note also that the spatial pattern in Fig. Etc) loses a 
periodicity of the original SkX, suggestive of a chaotic 
aspect of the melting dynamics. 

We finally compare the modes found in the SkX phase 
with those in the vortex-state nanodisks clarified in 
Rcfs. [26H29j . The twofold rotational modes and the 
breathing mode found in the SkX resemble, respectively, 
the twofold translational modes expressed by the Bcsscl 
functions with m = ±1 and the radial mode with m = 
in the vortex-state nanodisks. In Ref. (2||, Ivanov and 
Zaspel theoretically showed that degeneracy of the trans- 
lational modes with m = ±1 in the nanodisk is lifted un- 
der an applied magnetic field normal to the disk. We con- 
sider that a similar mechanism works in the SkX case for 
the doublet CW and CCW modes. There are also several 
differences. The modes in nanodisks are mainly governed 
by the long-range dipolar interaction, resulting in their 
salient aspect-ratio dependence. Note that their frequen- 
cies go to zero in the zero aspect-ratio limit. In contrast, 
the SkX and its dynamics considered here are governed 
by the nearest-neighbor spin interactions described in the 
Hamiltonian (1). The essential relevance of the DM in- 
teraction to the SkX is indicated by several experimental 
findings [IH such as its emergence only in chiral magnets, 
unique spin swirling directions of skyrmions, and consid- 
erably small size (10-100 nm) of skyrmions compared to 
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dipolar-forcc-induccd magnetic bubbles. Thus we expect 
negligible aspect-ratio dependence of the modes as well as 
weak influences of the dipolar interaction. Our study fo- 
cuses on thin films whose thickness is much smaller than 
the skyrmion diameter because a greater stability of the 
SkX in thinner films has been confirmed [13, [HI . In such 
a case, the system can be regarded as ferromagnetically 
stacked two-dimensional layers, which guarantees the va- 
lidity of our results based on a two-dimensional model. 

In summary, we have theoretically studied spin-wave 
excitations in the SkX phase of insulating fcrromagnets 
with DM interaction. We have found a couple of ro- 
tational modes with ~ 1 GHz frequencies for in-plane 
a.c. magnetic field. The rotations are in a CCW fashion 
for the lower-lying mode, while in a CW fashion for the 
higher-lying mode. These habits give rise to strong de- 
pendence of these spin-wave excitations on the circular 
polarization of the irradiating microwave. A breathing 
mode has been found for out-of-plane a.c. magnetic field. 
We have also observed the melting of the SkX under the 
irradiating LHP microwave. These findings will open a 
route to manipulation of skyrmion as a nano-scale spin 
texture using spin-wave resonances. 
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